Influence of a Random Telegraph Process on the Transport through a Point Contact 
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We describe the transport properties of a point contact under the influence of a classical two-level 
fluctuator. We employ a transfer matrix formalism allowing us to calculate arbitrary correlation 
functions of the stochastic process by mapping them on matrix products. The result is used to 
obtain the generating function of the full counting statistics of a classical point contact subject to 
a classical fluctuator, including extensions to a pair of two-level fluctuators as well as to a quantum 
point contact. We show that the noise in the quantum point contact is a sum of the (quantum) 
partitioning noise and the (classical) noise due to the two-level fluctuator. As a side result, we 
obtain the full counting statistics of a quantum point contact with time-dependent transmission 
probabilities. 
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I. INTRODUCTION 

The impact of mobile impurities on transport through 
a quantum conductor attracted a lot of attention soon 
after it was realized that the conductance of a dirty co- 
herent sample is sensitive to the position of a single im- 
purity; this discovery formed the basis for the explana- 
tion of nicker noise as it appears due to the presence 
of bistable mobile impurities. 1 -* 2 - A detailed characteriza- 
tion of charge transport through a quantum conductor 
is provided by the full counting statistics (FCS). 3 Dur- 
ing the past decade, this description has been applied to 
numerous systems 4 and a first attempt to describe the 
influence of a two-level fluctuator on the FCS of quan- 
tum transport has been given in Ref. 0. In the present 
paper, we calculate the full counting statistics of charge 
transport through a (classical or quantum) point contact 
coupled to a classical two-level fluctuator; we go beyond 
previous studies by considering the combined effects of 
the presence of one or many mobile impurities, as well 
as quantum-partitioning and the Fermi statistics on the 
FCS. Also, we reconsider carefully the case when parti- 
tioning is neglected and correct previous findings which 
are flawed when calculating the fourth or higher cumu- 
lants. 

The influence of a fluctuating (uncontrollable) envi- 
ronment on a (controllable) device is a generic problem 6 
and our work is related to other studies, e.g., the trans- 
port statistics through a quantum dot in the Coulomb- 
blockade regime*^ or the effect of a bistability on the 
transport through a quantum point contact^^. Another 
example is the dephasing of qubits due to a classical^ or 
quantum^ two-level system, or the studies on l//-noise 
originating from telegraph noise due to classica l 13 ! 14 or 
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two-level fluctuators. The two problems, 



full counting statistics of charge transport and dephasing 
of a quantum system (qubits) are related through the 
equivalence of fidelity^ and full counting statistics 3 -, as 
has been pointed out recently^ 

In our analysis below, we describe the time evolution 



of the two-level fluctuator by rate equations which can 
be solved explicitly. We then study the full counting 
statistics of a wire with a conductance depending on the 
state of the fluctuator. The fluctuator induces noise in 
the transport current through the wire which we evaluate 
using a mapping of correlation functions on matrix prod- 
ucts. Using this mapping, we are able to calculate the full 
counting statistics of a classical wire coupled to a two- 
level fluctuator. Furthermore, we discuss the situation of 
a second (independent) fluctuator and show that a non- 
linear interaction with the wire can lead to correlations 
in the noise even though the fluctuators evolve indepen- 
dent of each other. Finally, we apply our method to the 
case of a quantum wire which exhibits intrinsic partition- 
ing (shot) noise i 20 ' 21 We derive a formula which incorpo- 
rates both classical- (due to the two-level fluctuator) and 
quantum- (due to the point contact) noise. Thereby, we 
give an explicit expression for the full counting statistics 
of a quantum point contact whose transmission probabil- 
ities change with time. As two-level fluctuators seem to 
be a major obstacle for achieving solid state implemen- 
tations of qubits with long coherence times, being able 
to characterize the influence of a two-level fluctuator on 
transport through a point contact offers the possibility 
to learn about the fluctuating environment by measuring 
the full counting statistics through a nearby quantum 
point contact. The (partial) overlap of our results with 
previous worki^ii will be discussed below. 



II. SINGLE TWO-LEVEL FLUCTUATOR 

Consider a classical particle which can be trapped in an 
external potential at two positions denoted by x\ and x%, 
cf. Fig. [TJ The potential is characterized by the energies 
£1,2 associated with the two valleys, which are shifted 
by the amount A = E% — E\ with respect to each other, 
and the height of the barrier U. We assume the dynam- 
ics to be given by thermally activated hopping over the 
barrier. 22 The particle performs a random (Brownian) 
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FIG. 1: Sketch of the setup: A two- level system fluctuating 
incoherently between states £1,2 with rates 712 and 721 is 
coupled to a wire with a constriction. The conductance G 
of the wire changes according to the state of the fluctuator, 
inducing noise in the current flowing through the device. 

motion where the probabilities Pi j2 to be in either valley 
obey the rate equations 

A(i) = -72i-Pi(*) + 7i2-P 2 W, 

P 2 (t)= 721^1 (*) -712^0), (1) 

where the rate to escape from x± to X2 is given by 721 = 
7i->2 e" u ^ and the reverse process from X2 to x% is 
governed by the rate 712 = 72^1 oc e^ A ~ u ^^, with 1} 
the temperature. In equilibrium, the probabilities Pl\ 
are such as to obey the balance equation (IN21 = dN\2 
which equates the number of particles CL/V21 = ^lP^dt 
going from x\ to X2 during the time dt, with the number 
of particles dN 12 — JuP^dt going the opposite way. 
The equilibrium probabilities therefore satisfy the Gibbs 
weight 

PP/Pi q = e- A / § (2) 

and together with the probability conservation P^ q + 
P^ = 1, we obtain 

^ cq = Wr, P2 cq = 72i/r, (3) 

where we have introduced the total rate T = 712 + 721. 

Introducing the vector P(t) with components Pi^{t), 
the rate equation ([TJ can be written as P(i) = — hP(i), 
with the Fokker-Planck Hamiltonian^ given by 

h = ( 721 " 7l2 ^ . (4) 
\— 721 712 / 

Note that the Hamiltonian h is not Hermitian (left and 
right eigenvalues are not simply adjoint to each other). 
Nevertheless, its eigenvalues are real (and even positive). 
To make this point clear, we write the Hamiltonian h in 
a new basis h' = s _1 hs using the transformation matrix 
s = diag(-^/P° q , )' sucn that 

- {-r^pp 712 ) ■ (5) 



It is now visible that the matrix h' (and therefore also the 
matrix h) has real (the matrix is symmetric) and positive 
(the determinant and the trace of the matrix are positive) 
eigenvalues. The evolution conserves the probability as 
d t [Px{t) + P 2 {t)} = 0, as is evident from Eq. (TJ. This 
relation implies that (1,1) ■ P(i) = and it follows that 
(0| = (1,1) is a left eigenvector of the Hamiltonian h to 
the eigenvalue 0. To every left eigenvector there exists a 
corresponding right eigenvector with the same eigenvalue 
which we will denote by |0). The right eigenvector to 
the eigenvalue is given by the equilibrium distribution 
|0) = P cq , h 1 0) = 0. The second eigenvalue is given by T 
with the corresponding right [left] eigenvectors assuming 
the form |r) = (l,-l) T [(T\ = (P 2 eq , -P x eq )]; note that 
the eigenvectors are normalized such that (a\b) — 5 a b and 
J2 a \ a )( a \ = ^2 with a,b 6 {0,T}. Using the eigenbasis 
of h, it is possible to compute the evolution operator 
P(t > 0) = exp(-ht) = Ea e ~ at l a )H- The matrix ele- 
ment P m n(t) denotes the conditional probability for the 
particle to be transferred from state n to m in the time 
t. The evolution only depends on the time difference as 
the Hamiltonian h is time- independent. An explicit cal- 
culation yields the expression 

r[t) - y _ e -rt) P p + P p e -rt J W 

for the propagator during the time t. As the stochastic 
process is Markovian, the propagator ([5]) incorporates 
all the information needed in order to calculate general 
correlation functions of the stochastic process^ see also 
below. 



III. CLASSICAL WIRE 

We consider a classical wire coupled to the two-level 
fluctuator. We assume the two-level system to be a 
charge impurity which interacts with the wire, e.g., via 
Coulomb forces. The net effect of the charge impurity 
is to change the conductance of the wire depending 
on the position x(t) = x\$ of the two-level system. The 
wire is biased by a constant voltage V such that the cur- 
rent is determined by I± t 2 — VG± 2- The current in the 
wire jumps between I\ and J 2 in a random way given by 
the dynamics of the two-level fluctuator which we assume 
to be in thermal equilibrium, i.e., in the state |0). The 
fluctuations of the two-level system induce current noise. 
In our discussion, any kind of back-action of the wire on 
the two-level fluctuator is neglected. 

A. Correlation Functions — Mapping on Matrices 

With our focus on the full counting statistics, we 
are interested in obtaining the moments of the charge 
Q = L dtil(t\) transmitted through the point contact 
during the time t\ here, I(t±) denotes depending on 
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the state x(t\) — x\p of the two-level fluctuator at time 
t\. As I{t\) = VG(ti), we first concentrate on correla- 
tion functions of G(t), where statistical averages over the 
stochastic process (fTJ) will be denoted by (•). In quan- 
tum mechanics, it is well-known that correlation func- 
tions can be evaluated either in the operator or in the 
path-integral formalism. 25 Likewise, we have the choice 
to apply a stochastic path-integral approachi 2 ^ or to use 
the operator formalism. Here, we stay with the opera- 
tor formalism introduced in the previous section. Note 
that the propagator P mn (t 2 — ii), cf. Eq. ©, denotes the 
conditional probability (the transfer matrix) to find the 
system in state x m at time t 2l given that it resided in 
x n at time t\, P mn {t 2 — h) = (x(t 2 )=x m \x(ti)=x n )- 1 i.e., 
within a path-integral formulation, P(t) involves already 
an integration over all possible paths between t\ and t 2 . 

The average conductance (G(ti)) is readily calculated. 
For a system residing in a stationary state given by 
(x{t{)—x n ) — P° q , we obtain 

G(ti) = J2 G » <*(*o=*n> = G i p r + G * p 2 q - en 

n=l,2 

The calculation of the conductance correlator 
(G(t2)G(tx)) is more involved. We proceed slowly 
in order to motivate our general mapping between the 
calculation of correlation functions and the evaluation 
of matrix products of specific matrices. In order to 
calculate {G(t 2 )G(t 1 )), we assume first that t 2 > t\\ 
classical correlators are symmetric so that the opposite 
ordering of times reduces to the same quantity. Using 
the fact that the stochastic process is Markovian, 
we expand the correlation function (G{t 2 )G(tx)) = 
J2 m ,n G "i( x ( t 2)^x m \x(t 1 )^x n ) G n (x(t i )=x n ) ; 24 this 
expansion can be seen as a transfer-matrix expansion of 
the correlation function. We obtain the mapping for the 
correlator (t 2 > t\) 

(G{t 2 )G{h)) = G m P mn (t 2 - h)G n P^ 

mn, 

= G k SkmPml(t2-h)Gl6 ln P^ 
klmn 

= (0|Ge- h (* 2 - 4l )G|0) (8) 

with the diagonal matrix G rnn = G n 8 mn . Introducing 
the "interaction representation" 

Gi(f) = e ht Ge- ht (9) 

of the matrix G, the correlation function Eq. ([5]) can be 
further simplified to 

{G(t 2 )G{h)) = r<0|Gi(t 2 )Gi(ii)|0>, (10) 

where the time-ordering operator T has been included in 
order to relieve the restriction t 2 > t\. It is easy to see 
that the above derivation is not restricted to the second 
order correlation function, but can be applied in the same 



way to higher order correlation functions. We thus arrive 
at the mapping 

(G(t N ) ■ ■ ■ G(ti)) = T(0|Gi(tjv) • • • Gi(ti)|0), (11) 

where the left hand side is a correlation function for the 
classical stochastic process involving the two-level fluctu- 
ator and the right hand side is a matrix element involv- 
ing the matrices Gi(t„) and the vectors |0) = P cq and 
(OH (1,1). 

B. Full Counting Statistics 

We are now in the position to calculate the generating 
function 

x(A) = (e iA (12) 

for the zero-frequency current-correlation functions (mo- 
ments) of a classical point contact coupled to a two- 
level fluctuator. The moments are obtained as the Tay- 
lor coefficients (Q n ) — (— id\) n x\\=o- Alternatively, the 
stochastic process can be characterized by irreducible cu- 
mulants which are given by the expansion coefficient of 
the logarithm of the characteristic function 

((Q")) = (4r)" lo gx(A) a • (13) 

The characteristic function x(A) can be recast in the form 
X(A) = T(0|e av, £ di ' Gl (*')|0) (14) 

using the mapping Eq. (jlip . This equation can be 
rewritten using the well-known mapping between the 
"Schrodinger" and the "interaction" representation, 23 
e-( h+v )* = e- h *rexp[-/ *^'vi(t')] with vi(i) = 
e ht ve -ht Here, we apply the relation in the opposite 
direction to arrive at an expression without the awkward 
time-ordering, 

X (A) = (0|e ( - h+lAVG)t |0). (15) 

This formula was derived before by Bagrets and Nazarov 
using a stochastic path integral formulation of the prob- 
lem (the Fokker-Planck Hamiltonian h and the counting 
field A are denoted by L and A in their paper—). We 
believe that the present approach using transfer matri- 
ces and the mapping of the interaction picture onto the 
Schrodinger picture is more transparent. Note though, 
that the counting field A enters differently in their work 
compared to ours. Here, A couples to the classical current 
I whereas Bagrets and Nazarov discuss the transport of 
individual (quantum) particles such that A may only en- 
ter in the combination exp(iA) due to the quantization 
of charge. 

For further convenience, we subtract the average 
charge (Q) — V(G)t in order to obtain the reduced full 
counting statistics 

X(A) = X (A)e- lA(Q) = <0|e- fit |0), (16) 
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with the matrix h = h — i\V(G — (G)). Apart from the 
average charge (which is zero for x) , both log x and log x 
generate the same cumulants. The explicit calculation of 
the characteristic function of the full counting statistics 
X involves the eigenvalues 



1 + i\AgAP cq ± a/1 + 2i\AgAP^ - \ 2 (Ag)< 



of the matrix h, where we have introduced the difference 
in the equilibrium population AP cq = P^ — P± q and 
the difference in the (dimensionless) conductance Ag = 
V(G 2 -Gi)/T. 



which is different from the result 5 



((<3 4 )) 



*p?(p° 



peq )2 



V\G 1 - G 2 )H 

r 3 



The latter result is incorrect as it misses terms due to 
the implicit assumption in Ref. [H that the reduced con- 
ductance correlators {{G(t n ) ■ ■ ■ G(ti))) may only depend 
on the largest time difference t n — t\. Even though this 
hypothesis is correct for correlators up to n — 3, it fails 
for higher-order correlators. 



D. Short times 



C. Asymptotic long-time limit 

For long times Tt 3> 1, the matrix exponential (TTrll) is 
dominated by the eigenvalue of h with the smallest 
real part^ i 27 ' 28 We obtain an explicit expression for the 
generating function 



logX»0) = -ft-*- 



(17) 



All cumulants become linear in t, due to the fact that 
the autocorrelation time in the system is given by T _1 
and every state decays to the equilibrium state after this 
time. The fluctuations for Ft ^> 1 can be seen as a sum of 
independent stochastic processes and the cumulant gen- 
erating function logx becomes extensive in t. Interest- 
ingly, it is possible to obtain an explicit relation for the 
cumulants (n > 2) 



=n\V n (G 1 ~G 2 ) n T 1 - n t 

n-l 

xY, N ^A-V k+1 (Pi q ) k (P2 q Y 

k=l 



(18) 



with the Narayana numbers N n , m = ("J ( m ™ i) Tne 
cumulants in (|18j) grow factorially in magnitude with n 
and oscillate as a function of AP cq which are generic 
features of high-order cumulants^ The first couple of 
cumulants assume the form 



«Q 2 }>» = 



«q 3 »» -6Pr q p 2 cq (p 2 ° q -pr q 



r 



(19) 



,cqN V 3 (Gi - G 2 ) 3 t 



r 2 



Note that all the odd cumulants vanish if the process 
is symmetric 712 = 721, cf. Eq. (J3J. The cumulants in 
(HHJ) and HU) agree with the result in Ref. |5|. However, 
starting with the 4-th cumulant a discrepancy arises; e.g., 
for the 4-th cumulant we obtain 



«Q 4 » 



> 



24 peqpcq [( p0q _ peq )2 _ pOqpeqj 



VHGi-G 2 )H 

T 3 



For short times, t T^ 1 there is no evolution of the 
two-level system and we can set ht = 0. The full counting 
statistics reads 

^(A) = (0|e lAy(G - <G)) *|0) (20) 

= P 2 c V Ay(G2 ~ Gl)PlCqt + p ec l e ^^(<3i-G?2)P 2 e9 i 

with the cumulants (n > 2) given by 

((Q n ))«=V n (G 1 -G 2 ) n t n (21) 

;E n - 1 , k - 1 (-i) k+1 (p? l ) k (i? l ) n - k , 



n-l 



k=l 



where the Eulerian numbers E n ,m are defined through 

E n, m = Er=o(- 1 ) fe (™fc 1 )( m + 1 - A: )"'^ Intne SnOTt time 

limit, the cumulants grow like {(Q n )) oc t n , i.e., higher 
order cumulants are suppressed at short times. The first 
couple of cumulants are explicitly given by 



{(Q 2 ))»-^r^ oq ^(Gi-G2) 2 t 2 , 

«Q 3 »» - Pr q P2° q (^2 q - P° 1 q )V 3 (G 1 - G 2 ft 3 



(22) 



E. Arbitrary times 

Expanding the matrix exponential in Eq. (|16[) in its 
eigenbasis, the generator for the full counting statistics 
reads 



-h-t 



h-e 



-h+t 



h + — h- 



(23) 



for arbitrary times. This result has been first derived in 
Ref. [ll] in the context of dephasing of a qubit due to the 
interaction with a classical two-level fluctuator, where A 
denotes the interaction of the qubit with the fluctuator. 
Here, we are interested in the transport properties of a 
point contact characterized by cumulants which are given 
by the Taylor expansion of logx around A = 0; the rela- 
tion between these two problems is a consequence of the 
generic equivalence between full counting statistics and 
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fidelity, see Ref. lljj The first couple of cumulants are 
given by 



«Q 2 }} = 2P eq P 2 ° q 



V 2 (Gi-G 2 ) 2 {(Tt-l) + e 

r 2 



-rti 



((Q 3 )) = 6P° q PP(Pp - Pr q ) (24) 
^ 3 (Gi - G 2 ) 3 [(n - 2) + (R + 2)e- rt ] 

x fs ■ 



F. Symmetric fluctuator 

A special situation is given when the two-level fluctu- 
ator is symmetric, A = 0, i.e., = P^ q — 1/2. Then 
the characteristic function assumes the simple form 

Tf 



logX»(A) = — [y/i - A 2 F 2 (C - G 2 ) 2 /r 2 - lj (25) 

for long times. In the short time limit, the generating 
function 



X<(A) = cos[AV(Gx-G 2 )V2] 



(26) 



becomes periodic. In both cases, due to the symmetry 
of the states x\ and x 2 , only the even cumulants are 
nonvanishing. 



IV. A PAIR OF TWO-LEVEL FLUCTUATORS 



Needless to say, the mapping of Sec. IIII Al is not re- 
stricted to a single two-level fluctuator. It can be gen- 
eralized to an arbitrary number of states whose dynam- 
ics is governed by classical rate equations described by 
a Fokker-Planck Hamiltonian h. To illustrate this con- 
cept, the example of two classical, uncorrelated two-level 
fiuctuators coupled to a wire is discussed in the follow- 
ing. We restrict ourselves to the case where the dy- 
namics of the two-level systems is completely indepen- 
dent of each other such that h = h" + h* 3 ; here and 
in the following, we denote quantities involving only 
the first (second) fluctuator with a superscript a(/3), 
e.g., h" = h" ® flr. Written explicitly in the basis 
{|1}" ® |1}' 3 ,|2)" ® |l)Ml) a ® |2)< 9 ,|2)" g> |2)' 9 }, the 
Fokker-Planck Hamiltonian reads 



\ 



-721 

-721 





7l2 



-712 



"721 



-7f 2 



721 



712 





-7? a 



-721 



-712 

7i°2+7iV 



(27) 



One is tempted to think that the independent dynam- 
ics of the two subsystems would lead to independent 
statistics, such that the characteristic function of the 
full counting statistics is given by the product of the 
individual characteristic functions. Indeed, this is the 
generic case for two-level fiuctuators coupling to a qubit, 
where the combined effect leads to 1/f noise i 11 ^ 3 ' 15 " — 



However, here, this argument is only valid when the in- 
teraction with the wire is "linear" such that the effects 
of the individual subsystems simply add up, in formula 
G = G" + G@. Having the model of the Fig. Q]in mind, 
this assumption is incorrect as a (quantum) point contact 
does not react linearly on changes in the gate potential 
and, therefore, the noise of individual fiuctuators does 
not simply add up. In the following, we first treat the 
(simple) case of linear interaction and then comment on 
the correlation which arises in the general case by apply- 
ing perturbation theory in the nonlinearity. 

Introducing the reference conductance Go = Gq + Gq 
as well as the induced changes AG X = Gf — G x due to 
the fluctuator x = a, /3, the conductance matrix for linear 
interaction is given by 

G = G 1 4 + diag(0, AG", AC 3 , AG" + AC 3 ); (28) 

the increase of conductance in the case when both fiuctu- 
ators are in state x 2 is the sum of the respective increases 
when only one of the fiuctuators is in state x 2 , that is to 
say, the effects of the two fiuctuators simply add up. In 
this case, the characteristic function assumes the form 
x(A) = X a (A)x /3 (A) and the cumulants ((Q n )) become a 
sum of cumulants generated by the two individual sub- 
systems. 

In the case of a general (diagonal) matrix G, the so- 
lution of the problem involves the determination of the 
roots of a polynomial of fourth degree and the charac- 
teristic function does not separate, even though the time 
evolution of the two fiuctuators is completely indepen- 
dent of each other. To be more explicit, we want to show 
how a small perturbation AG <C AG" + AG^ in G44 
destroying the linearity (additivity) leads to correlations 
which can be arbitrary large for long times. We define 
x co "(X) = x(A)/x Q (A)x /3 (A) as the part of the charac- 
teristic function which describes the correlation between 
the action of the individual subsystems. In the long- 
time limit, to first order in AG, we can apply standard 
perturbation theory to find the correction to the lowest 
eigenvalue of Eq. (|27[) . The cumulant generating function 
for the correlation is given by 

logxrxA) = 

-n- (T x + i\VAG x )PP' x - h x 
iXVAGt Y[ -■ (29) 



x—a,0 



h x , - h x 



The average transmitted charge changes according to 



A(Q) = P 2 cq '"P 2 eq ' /3 yAGt, 



(30) 



with p^pp'P the probability to be in the state |2}" ® 
|2}^ and VAG the change in the current. Less trivial, 
the correlation contribution to the noise 



A((Q 2 }} = WA(Q) 



peq,a 
pa 



AG" 



PI 



cq,/3 



re 



-AG? 



(31) 



depends both on AG" and AC 3 . 
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V. QUANTUM WIRE 

Considering a quantum rather than a classical wire, 
additional noise appears due to the probabilistic nature 
of the charge transport (due to partitioning) even in the 
absence of a fluctuating environment. Given a quantum 
wire with N channels characterized by their transmission 
eigenvalues T 7 , 7 = l,..,,N, and biased by a voltage po- 
tential V, the characteristic function of the full counting 
statistics is given by£ 



logX q (A) = |^ ^ log [1 + (e** - 1)T 7 ], 



(32) 



with q the charge of the electron; this result is valid in the 
asymptotic limit qVt/h 3> 1, for low-temperatures i? <C 
ql^, and with the proviso that the energy-dependence of 
the transmission eigenvalues is negligible in the energy 
interval set by the voltage. The quantum nature of the 
fcrmions leads to the noise 



«Q 2 » q = q(Q) q 



E 7 r 7 (i-r 7 ) 



(33) 



which disappears provided that all the channels are either 
closed T 1 — or completely open T 7 = 1. Note that the 
noise in Eq. (|33|) is sub-Poissonian, i.e., the Fano factor 
F = ((Q 2 )) q /q(Q) q is less than 1. 

Here, we are interested in the case where the quan- 
tum wire is capacitively coupled to a two-level fluctuator 
such that the transmission eigenvalues T 7 change over 
time; note that we neglect a possible energy dependence 
of the transmission eigenvalues, which corresponds to the 
fact that we assume that the scattering center does not 
produce any time delay due to the scattering event. The 
characteristic function logx q (A) = det Q is given by the 
determinant of the matrixsl 



Qfc 7 ,fc' 7 ' = {4> kl {t)\e lXqQt \4> klll {t)) 



(34) 



with the counting operator Q t = Jrdx \x)(x\ integrated 
over the interval I = [0, v F t] and 



(j) kl (x,t) 



' e ik(x-v v t) +Ti (t + x/v F )e- lk{x+v ^ x < 
r^t- x/v F )e lk{x - VFt) 



x > 
(35) 

the single-particle solution of the time-dependent 
Schrodinger equation involving a scattering center at 
x = with time-dependent transmission [reflection] am- 
plitude t 7 (£) [r 7 (i)]; note that we have suppressed the 
transverse part of the wave function belonging to the 
channel index 7. Equation (I35p is valid in the linear- 
spectrum approximation, where v F is the Fermi velocity 
and < k < k F 2^ As the matrix Q is block-diagonal in 7 
and constitutes a Tocplitz matrix with respect to the in- 
dex fc, its determinant can be shown (using the technique 
discussed in Ref. Hl|) to have the form 

logX q (A) = g fdt> 5>g[l + (e^-l)T 7 (i')] (36) 



for long times qVt/H 1; the expression Eq. (l36l) reduces 
to Eq. (|32j) when the transmission probability does not 
change in time. We now add the classical two-level fluc- 
tuator to the system: Depending on the state x — x\^ of 
the nearby two-level fluctuator, the quantum wire is de- 
scribed by one of the two sets of transmission eigenvalues 
Ti 2 - The total generating function x qtl (A) is an average 
over the individual contributions 



X qtl (A) = (x q (A)) 



= (JqV/lirh) £ dt' £ 7 log[l+(e < ^-l)T 7 (i')] 



>, (37) 



where the average (•) is over the stochastic process of 
the two-level fluctuator. Equation (|37l) can be calculated 
explicitly with the method outlined in Sec. lIII Al Indeed, 
Eq. (TT2"j) goes over to Eq. (|3"T)l via replacing iXI(t) with 
(qF/27rft)£ 7 log[l-Ke 4 « A -l)T 7 (f)]. In this mapping, the 
conductances G n are changed to (q/2nihX) ^2 log[l + 
{e tqX — 1)T 7 ] with T 7 the transparency of the channel 7 
when the two-level fluctuator is in the state n. Inserting 



iXAg — > /j 



53 log 



1 + (e l « A - 1)T 2 7 
1 + (e l « A - I)!? 



(38) 



into Eq. (fTT)) and using Eq. (ITBl . we obtain 
logX^>(A)=5>riogx;?(A) (39) 



n,7 

rt 

y 



1 + n AP cq - Vl + 2^AP°q + /i 2 



here, log X «(A) = (qVt/2irh) E 7 log[l + (e^ A - 1)T 7 ] is 
the characteristic function for the FCS, Eq. (|3"2j) . depen- 
dent on the system's state (n = 1,2) and 7 runs over the 
channel index. The expression (|39[) coincides with the re- 
sult obtained by Jordan and Sukhorukov who considered 
the case of rare transitions (HT <c qV), see Ref. [9J. In 
their work, Jordan and Sukhorukov describe the trans- 
port of a conserved charge in a classical bistable system 
where large charge fluctuations drive transitions between 
stable states with different transport characteristics. The 
generic fluctuations present in the system's stable states 
then combine with the fluctuations of the bistable charge 
to generate the transport statistics of the bistable sys- 
tem. While describing a very general situation, the spe- 
cific analysis in Ref. [9| is limited to those cases where the 
charge-switching events are rare on the time scale of the 
underlying fluctuations in the stable states. The transla- 
tion to our model may not be obvious from the start and 
is done by choosing the process of charge partitioning 
for generating the underlying fluctuations (with different 
transmission rates T 7 , n = 1, 2 for the two stable states), 
i.e., the generators logx2,(A) correspond to the long time 
generators H n t in Ref. [9|. The non-linearity leading to 
switching between stable states generates the transition 
probabilities Ti^; the latter arc determined by an instan- 
ton trajectory and have to be small in the case of Ref. 0. 
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qV/hr 



FIG. 2: Plot of the Fano factor F as a function of the bias 
voltage V for a single mode wire with Pj q = = 1/2, 
Ti = 1, and T 2 = 0, 0.1, 0.3, and 0.5 . The Fano factor starts 
of at a value T 2 (l - T 2 )/(l + T 2 ) < 1 for small voltages. In the 
opposite regime, it is given by [(1 - T 2 ) 2 /(l + T 2 )]qV/2nhr 
which becomes superpoissonian for V large enough. 



In our model these rates are given as the basic input pa- 
rameters 7ij defining our two-level fiuctuator. The con- 
dition of rare charge switching events corresponds to the 
requirement that fir <C qV. In this limit Eq. (|39p reduces 
to the result of the classical point contact, Eq. (|17p . with 
the conductances given by 



g 

2nh 



(40) 



We observe that the statistics is dominated by the fluc- 
tuations of the impurity and the results derived in the 
previous sections remain valid in the case of a quantum 
wire when the classical conductance is replaced by the 
Landauer formula f4"0")) . 

However, our result Eq. (f39|) is also valid in the op- 
posite regime fiF 3> qV, i.e., when the two-level fiuctu- 
ator noise acts on a timescale which is fast compared to 
the partitioning noise. Then only the first term in (|39[) 
contributes and the cumulant generating function is the 
average of the expressions (|32[) for the quantum point 
contact, to be taken over the positions x\^ 2 with weights 
given by the probabilities P^\- Having access to both 
regimes, it is possible to study the crossover from clas- 
sical noise (due to the two-level fiuctuator) to quantum- 
partitioning noise (due to the quantum point contact). 
To this end, we calculate the long-time asymptotics of 
the first two moments of (|3"T|) , 



(Q)f = V(P^G 1 + PpG 2 )t 
for the average charge and 



(41) 



V 2 {G 1 -G 2 )H 



+ 2P 1 Cq P; 



eq 
2 " 



(42) 



for the noise. The noise is simply given by the sum of 
the quantum partitioning noise (first term) and the noise 
due to the dynamics of the impurity (second term) . Note 
the crossover of the noise (l4"2l from sub-Poissonian F < 1 
for a fast fiuctuator with fir ^> qV [the first term in (|42)) 
dominates] to super-Poissonian F > 1 when the fiuctua- 
tor is slow fir <C qV [the second term in (|42"T) dominates] , 
provided that G\ ^ G 2 , see Fig. [2] 



VI. EXPERIMENTAL TEST 

It is difficult to experimentally confirm the crossover 
from sub- to super-Poissonian noise as described in the 
previous section as the quantum-partitioning noise is typ- 
ically small and thus the noise of the classical two-level 
fiuctuator will dominate. As a promising setup, we envi- 
sion coupling a quantum point contact in GaAs/AlGaAs 
to a double dot, e.g., in an InAs nanowire, which serves as 
tunable two-level fiuctuator. Such a system was studied 
recently, see Ref. |H. 

To observe the crossover, both the noise Sti(0) from 
the classical two-level fiuctuator and the quantum- 
partitioning noise 5 qp (0) (at zero frequency) have to 
dominate over the thermal Nyquist- Johnson noise which 
is given by 



Snj(0) 



2g 2 k B T 



1CT 



; A 2 s 



(43) 



at 50 mK. The quantum-partitioning noise [first term in 
Eq. (|4"2")) ] can be estimated as 



«M0) 



q 3 V 
8^fi 



V"[mV]10~ 27 A 2 s 



(44) 



when the system is tuned in the middle of a conductance 
step with T w 1/2. Note that for bias voltages V > 
0.1 mV the quantum-partitioning noise is larger than the 
thermal noise floor. 

Tunneling of an electron between the two quantum 
dots with a rate T leads to a change in the conduc- 
tance. In Ref. l33l. this change was of the order of 0.1 q 2 /h. 
However, in order to being able to observe the crossover 
the capacitive coupling of the double dot to the quan- 
tum point contact should be reduced to a level such that 
G\—G 2 ~ 0.001 q 2 /h. This provides us with the estimate 
(using PP = PP = 1/2) 



Sti(0) = 



V 2 (G 1 -G 2 ) 2 _ y 2 [mV 2 ] 

r ~ r[MHz] 



10" 



; A 2 s (45) 



for the zero-frequency noise due to the two-level fiuc- 
tuator [second term in Eq. (j42|)]. At a bias voltage 
V ~ 1 mV with a rate r ~ 100 MHz, the noise due to the 
fiuctuator is given by Sti(0) ~ 10~ 28 A 2 s dominated by 
the quantum-partitioning noise with Sti(0) ~ 10~ 27 A 2 s. 
Note that in experiments rates of the order of 10—100 kHz 
have been observed j 33 ' 34 We expect that rates in the 
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100 MHz regime to be realistic due to the exponential 
dependence of the tunneling rate on the potential bar- 
rier. Increasing the bias voltage V to values of a few mV 
the noise due to the classical two-level fluctuator starts 
to dominate and the crossover as depicted in Fig. [2] can 
be observed. 



VII. CONCLUSION 

We have determined the influence of a thermally driven 
two-level fluctuator on a point contact through the cal- 
culation of the full counting statistics of transported 
charge. Both, classical and quantum point contacts have 
been considered and extensions to multiple fluctuators 
have been discussed. In our analysis, we have made 
use of a mapping between correlation functions of classi- 
cal stochastic processes and simple time-ordered matrix 



products. For the case of a quantum point contact, we 
have shown that the partitioning noise and the noise due 
the two-level fluctuator add up and the noise crosses over 
from sub- to super-Poissonian depending on the applied 
voltage bias. To extend the present formalism to the case 
of current correlators at finite frequencies, which provides 
additional insights into the dynamics of the two- level fluc- 
tuator, is an interesting problem for future studies. 
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